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quantify the quantum gravitational back-reaction on inflation. Although our 
dynamical assumption may not prove to be completely correct, it does have 
5^ ! the right time dependence, it can naturally produce primordial perturbations 

of the right strength, and it illustrates how a rigorous application of the 
leading logarithm approximation might work in quantum gravity. It also 
serves as a partial test of the "null hypothesis" that there are no significant 
effects from infrared gravitons. 
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1 Introduction 

There are profound differences between quantum gravity with a positive cos- 
mological constant and its flat space version . Many of these differences are 
associated with the phenomenon of infrared logarithms. These are factors 
of the number of e-foldings since the onset of inflation that show up in the 
Green's functions of theories which incorporate either gravitons or massless, 
minimally coupled scalars. Their physical origin is the continual creation 
of long wavelength gravitons and scalars during inflation, which engenders a 
slow growth in the corresponding average field strengths. Infrared logarithms 
have been seen in pure quantum gravity [T] , in gravity + fermions [2] , in full 
scalar-driven inflation [3111], in the scalar sector of scalar-driven inflation [S], 
in ip^ theory [6j , in scalar QED [7] and in Yukawa theory [8], [9] . 

Infrared logarithms introduce a fascinating, secular element into the usual 
static results of flat space quantum field theory. For example, whereas the 
virtual gravitons of flat space have no net effect on the propagation of a single 
fermion, the spin-spin coupling with the sea of infrared gravitons produced 
during inflation makes such a fermion behave as if its field strength were 



being amplified by a factor of l/JZ2{t) with [2]: 

Z2{t) = 1 -^GH'J^]+0{G') . (1) 

Here G is the Newtonian gravitational constant, H the inflationary Hubble 
parameter and Oon represents the scale factor at the onset of inflation. The 
dimensionless product GH^ < 10~^^ is bounded by the current upper limit 
on the tensor-to-scalar ratio [TDj. Despite the minuscule coupling constant, 
the continued growth of the scale factor a{t) over a prolonged period of in- 
flation must eventually result in this 1-loop correction becoming order unity. 
Because higher loop corrections also become order unity at about the same 
time it is clear that perturbation theory breaks down and that one must 
employ a nonperturbative technique to evolve further. 

The eventual breakdown of perturbation theory evident in expression ([1]) 
is a general feature of quantum field theories that exhibit infrared logarithms. 
Starobinskii has developed a simple stochastic formalism [11] which has been 
proved to reproduce the leading infrared logarithms of scalar potential models 
at arbitrary loop order [12]. When the scalar potential is bounded below this 
technique can even be used to evolve past the breakdown of perturbation 
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theory to asymptotically late times P^. Starobinskii's technique has recently 
been generalized to scalar models which involve fields that do not produce 
infrared logarithms such as fermions [9] and photons [H]. However, it has 
not yet been extended to the derivative couplings of quantum gravity. 

The purpose of this paper is to explore quantum gravity at leading loga- 
rithm order. Q| We have two results, one rigorous and the other qualitative. 
The rigorous result is that "leading log" quantum gravity is much better 
behaved than full quantum gravity. The only counterterm that contributes 
at leading logarithm order is the 1-loop renormalization of the cosmological 
constant. This is derived in Section 2. 

Our qualitative result consists of a simplifying dynamical assumption, 
the Effective Scale Factor Approximation, under which the leading infrared 
logarithms of quantum gravity can be computed and summed to give non- 
perturbative results. Although this assumption may not be entirely right, its 
qualitative features are consistent with all that is currently known about in- 
frared logarithms, and the approximation allows us to explore how a rigorous 
application of Starobinskii's formalism might work in pure quantum gravity. 

Both results stand on their own, however, it is interesting to view them 
in the context of our long-held suspicion that infiation might be driven by a 
positive bare cosmological constant without any scalar infiaton [T71 [1] . An 
attractive feature of this idea is that it dispenses with the usual problem of 
explaining why the cosmological constant is so small by assuming that A is 
actually GUT-scale. That would start infiation in the early universe without 
the need to assume an unnaturally smooth initial condition. Infiation would 
be stopped in this model by the back-reaction from long wavelength, virtual 
gravitons which are continually ripped out of the vacuum by the accelerated 
expansion of infiation. The kinetic energy density of these gravitons is ex- 
actly diluted by the expansion of the 3-volume to produce a constant energy 
density. However, the next order effect, due to the gravitational interaction 
energy between gravitons, must slow infiation because gravity is attractive. 
This next order effect grows without bound, so it must eventually stop infla- 
tion if nothing else supervenes flrst. 

One can understand the growth of back-reaction from the fact that the 
Newtonian potential is —GM/R, with the total kinetic energy of infrared 
gravitons growing like M ~ a^{t) and the radius growing like R ~ a{t). Hence 
the universe must eventually fall inside its own Schwarzschild radius! That 



^For some interesting alternate approaches, see |15[I16). 



estimate (of a^(t) growth for the energy density of interaction) assumes the 
newly created gravitons are instantly in contact with one another, whereas 
the self-interaction actually requires time to build up. When this causality de- 
lay is taken into account the interaction energy density grows like —GH^ ln(a) 
at lowest order [18]. This slower growth means that inflation can only be 
checked after an enormous number of e-foldings: ln(a) ~ 1/{GH'^) ~ 10^^. 
There seems to be no problem with such a large number of e-foldings; indeed, 
it can be viewed as a way of exploiting the weakness of the gravitational in- 
teraction to make inflation last a long time. We will see that it is also a way 
of keeping spatial variation small even over volumes as large as the currently 
observable universe. 

Although we emphasize that the results of the present work stand on their 
own, they can be viewed as a partial check on the "null hypothesis" that 
infrared gravitons make no significant corrections to inflationary quantum 
gravity [TH [20]. Further, if one assumes that the spatially homogeneous, 
leading infrared logarithms stop inflation, then the spatially inhomogeneous 
corrections at the next order are suppressed by one factor of GH^, which is 
the observed strength of the spectrum of primordial density perturbations. 
Hence scalar density perturbations of the right strength seem to be possible 
without a fundamental scalar. 

Section 3 motivates and defines the Effective Scale Factor Approximation. 
In Section 4 we compute the effective scale factor, both perturbatively using 
dimensional regularization at 1- and 2-loop orders and, implicitly, at arbitrary 
order. Section 5 sees this result used to evaluate an invariant observable 
which has been proposed to quantify the quantum gravitational back-reaction 
on inflation [21]. Our conclusions comprise Section 6. 

Before closing this Introduction we should clarify the distinction between 
"full quantum gravity" and "quantum gravity at leading logarithm order". 
Because general relativity is not perturbatively renormalizable, it cannot pro- 
vide a complete theory of quantum gravity on the perturbative level. There 
must either be some different model, possibly not even based upon a metric 
field, or else quantum gravity computations are intrinsically nonperturba- 
tive. However, failing to determine everything is not quite the same thing 
as failing to determine anything. It is perfectly valid to employ perturbative 
quantum general relativity as an effective field theory to study phenomena 
which are driven by infrared gravitons. Further, the results so obtained can- 
not be changed by the still unknown, ultraviolet completion of the theory. 



So we will many times derive results which take the form, 

(finite number) x ln(a) + (divergent constant) , (2) 



and retain the infrared logarithm while ignoring the divergent constant. This 
is correct because the infrared logarithm derives exclusively from long wave- 
length virtual gravitons, which must be reliably described by quantum gen- 
eral relativity. In contrast, the divergent constant originates in the ultraviolet 
sector, which cannot be correct, at least not on the perturbative level. Be- 
cause no divergences are observed in nature we know that the ultraviolet 
completion of quantum gravity is somehow avoiding them. As long as the 
result is a finite constant, the infrared logarithm must eventually dominate 
at late times. 

This use of a fiawed (or misunderstood) formalism as an effective quan- 
tum field theory to infer valid results from the infrared has a long and distin- 
guished history. The oldest example is the solution of the infrared problem 
in quantum electrodynamics by Bloch and Nordsieck [22], long before that 
theory's renormalizability was suspected. Weinberg [23] was able to achieve a 
similar resolution for quantum gravity with zero cosmological constant. The 
same principle was at work in the Fermi theory computation of the long range 
force due to loops of massless neutrinos by Feinberg and Sucher [23] . Matter 
which is not supersymmetric generates nonrenormalizable corrections to the 
graviton propagator at one loop, but this did not prevent the computation of 
photon, massless neutrino and massless, conformally coupled scalar loop cor- 
rections to the long range gravitational force [25]. More recently, Donoghue 
[26j has touched off a minor industry [27] by applying the principles of low 
energy effective field theory to compute graviton corrections to the long range 
gravitational force. Our analysis exploits the power of low energy effective 
field theory in the same way, differing from the previous examples only in 
the detail that our background geometry is locally de Sitter rather than flat. 

2 Quantum Gravity at Leading Log Order 

The purpose of this section is to describe the leading logarithm approxima- 
tion for inflationary quantum gravity and to demonstrate that the theory in 
this limit is vastly better behaved than full quantum gravity. We begin by 
reviewing the free field expansion of quantum gravity on a locally de Sitter 
background. We then take note of some key facts about infrared logarithms. 



Finally, we show that only the 1-loop renormalization of the cosmological 
constant can contribute at leading logarithm order. 

2.1 The Free Field Expansion 

Because we employ dimensional regularization it is necessary to work in D 
spacetime dimensions. With a cosmological constant A, the gravitational 
equations of motion are: g 



R^,.--g^uR+-{D-2)Ag^, = Q 



(3) 



The free field expansion of perturbative quantum gravity has much the same 
structure on a locally de Sitter background as on flat space. Consider the 
case of synchronous gauge for which the invariant element takes the form: 



ds^ = —dt^ + c?{t) exp K/i(t, x) . . dx'dx^ 



»j 



a{t) = e 



Ht 



(4) 



Here, the exponential of the graviton field /ijj(t, x) has its usual meaning: 

1 



exp 



nh 



V 2 



ik'^kj 



(5) 



where k^ = 16nG is the loop counting parameter. 

Just as in fiat space we can decompose the graviton field into irreducible 
representations of the rotation group: 
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h- = /if-' 4-Bh +dh 



S^J-{D-l 



didj 



h^ + 



1 



D-2 



didj 



V2 



h , (6) 



where the usual conditions distinguish the transverse-traceless and transverse 
components: 

d.i^;" = hr = uj = . (7) 

Like fiat space, the linearized equations of motion imply: 

hj{t,x) = h^{t,x) = h{t,x) = , 



(8) 



^Hellenic indices take on spacetime values while Latin indices take on the {D — 1) 
space values. Our metric tensor has spacelike signature and our curvature tensor equals 



R"pt.u = ^".p,^. + r"^p r"^^ - (a* ^ ^)- The Hubble parameter is H : 



S^A. 



while the hnearized solution for the transverse-traceless components can be 
expanded in spatial plane waves: 

X^,it, x) ^ v^l A-^ eik-H) Y: { uit, k) e'%,ik, A) aik, A) + (cc) } . 

(9) 
In dn]), u(t,k) are the mode functions, eij{k,\) is the polarization tensor, 
a{k,X) the annihilation operator, and (cc) stands for complex conjugation. 
The |(Z) — 3)D creation and annihilation operators and transverse-traceless 
polarization tensors obey the same relations as in flat space 



[a{k,X),a\k',X')] = {27rf-'6xx'S''-\k-k') , (10) 

J2e,,ik,\)eUk,\) = n,^kik) n,)^(k) - ^^n,,(k) U,,(k) , (11) 

A 

eu{k,X) = kieij{k,X) = , ei,(^. A) e^(^. A') = 6xy , (12) 

where Ilij{k) is the transverse projection operator: 

iijjl^fcj = Oij fcj Kj , (,l"^j 

and where parenthesized indices are symmetrized. 

The main difference from fiat space is that the mode functions u{t, k) go 
from simple exponentials to Hankel functions: u 

„,a)^>;j|.-<v.^«,(_L). (14) 

The restriction to k = \\k\\ > H in iQ is made to avoid an infrared singularity 
in the free propagator [3U] . The physical reason for this singularity is that no 
causal process would allow a local observer to prepare coherent Bunch-Davies 
vacuum over an infinite spatial section. Sensible physics can be regained 
either by employing an initial state for which the super-horizon modes are 
less strongly correlated [31] , or else by worlcing on a compact spatial manifold 
such as T^~^ for which there are initially no super-horizon modes [32], and 
then making the integral approximation to the mode sums. In both cases 
modes with k < H are effectively absent. 



^The system is assumed to be released in free Bunch-Davies vacuum at t — 0. 



The parallel with flat space persists at higher orders. The fundamental 
quantity of the free field expansion is the linearized mode sum (^. All 
components of hij can be expressed in terms of Xij- The qqq and goi constraint 
equations determine hj , h^ and h as expansions which start at second order: 

hj ^h^ ^h^ Kx^ + K^X' + . . . (15) 

The dynamical gij equations define a related expansion for hj^ which starts 
at first order: 

hj^ ~ X^, + /^X' + /^'x' + • ■ ■ (16) 

These expansions will generally involve integrations against retarded Green's 
functions - which can themselves be expressed as the commutator of two free 
fields: 



ijGke 



{x;x') = ie{t-t') Xij{x),Xke{x') . (17) 



A nice diagrammatic representation for the free field expansions has been 
given recently by Musso [5^ . 

The expectation value of any operator can be computed using the free 
field expansion. One first expands the operator in powers of the graviton 
field hij. Each term in this series gives its own expansion in powers of the 
free field Xij- O^is then takes the expectation value of the resulting sums 
and integrals of products of free fields, for which the usual reductions of free 
field correlators apply. In particular, the expectation value of N of the Xij's 
is zero for N odd, while for A^ even it is the sum of the (A^ — 1)!! distinct 
permutations of 2-point correlators. 

2.2 Facts about Infrared Logarithms 

We have already seen the crucial role of infrared logarithms for perturbative 
quantum gravity on a locally de Sitter background. These infrared logarithms 
derive from explicit factors of ln(a) □ in free field correlators [M], and from 
integrating the retarded Green's functions (times these correlators) that arise 
in the free field expansion [5SII5]. Based on this understanding one can give 
a simple rule - valid for any interaction in any theory - for counting the 
maximum number of infrared logarithms that can arise for each extra power 



*To simplify the notation we have normahzed the initial value of the scale factor to be 
= 1. 



of the coupling constant. Consider an interaction with A^ undifferentiated 
massless, minimally coupled scalars or gravitons: 



A£ ~ (coupling constant) x {(p,hij) x (anything 



(18) 



Here, "anything" can involve differentiated if or /ijj's, or other fields alto- 
gether. Then, for each additional factor of the square of the coupling constant 
there can be at most A^ additional infrared logarithms [H]. A few examples 
from scalar models on a non-dynamical de Sitter background illustrate the 

rule ^nm- 



->^a- 



AC 



AC 

4! 



AC 



-fipipip a 



D 



Aln2(a) 
e^ln(a) 



(19) 

(20) 
(21) 



The basic interaction of the bare quantum gravitational Lagrangian has 
long been known for de Sitter background [36] and it gives the result [H] : 

(22) 



AC ~ Khdhdha^-^ =^ GH^ln{a) . 
The counting is the same for gravity + fermions [2]: 

AC ~ Khiidi)a^-^ =^ GH^\n{a) . 
Hence, the general form of the fermion field strength ([1]) is: 



(23) 



Z2it) 



1 + Y,{GH' 



Q,o 



In (a) 



+ Q,i 



ln(a) +... + Q,^_iln(a)[ . (24) 



The constants C£,k are pure numbers which are assumed to be of order one. 
The term in flMl) involving [GH"^ ln(a)]^ is the leading logarithm contribution 
at i loop order; the other terms are suhdominant logarithms. 

Quantum gravitational perturbation theory breaks down when In (a) ~ 
l/GH"^, at which point the leading infrared logarithms at each loop order 
contribute numbers of order one. In contrast, the subleading logarithms are 
all suppressed by at least one factor of the very small number GH'^ < 10^^^. 
So it makes sense to retain only the leading infrared logarithms: 



Z2{t) 



l + Y,ceJGHHn{a) 



(25) 



This is known as the leading logarithm approximation. 

It is important to note that, the operator under study can affect the form 
of its leading logarithm expansion. For example, we saw from expression 
flT9|) that one gets at most two infrared logarithms for each extra coupling 
constant from a Xip^a^ interaction. That rule only gives the maximum, and 
it does not specify the number of infrared logarithms in the lowest order 
result. Explicit computation shows that the leading logarithm expansion of 
the expectation value of 2n coincident fields is p^: U 



((fi|v?'"(x)|fi)) = (2n-l)!! 



H"^ ln(a) 



47r2 
35n^ + 170^3 + 225n2 + 74n 



n^+n A 2/ N 
ma 



2 367r2 
+ 



280 



^ In^fa) 



367r2 



(26) 



In this case the bound of In^(a) for each extra A is saturated for every n, but 
the In (a) dependence of the order A° result depends upon n. In contrast, the 
expectation value of the kinetic term takes the form [6] : 

{{n\g^''d,^d,^\n)) = -^|l + #Aln(a) + 0(A2ln3(a))| , (27) 

where # is a divergent constant jj Note that the derivatives in the operator 
under study not only preclude there being any infrared logarithms at zeroth 
order, they also cause loop corrections to be one fewer than the maximum 
given by the rule flT^ . So instead of the order A correction possessing two 
infrared logarithms, it has only one; instead of the order A^ correction having 
four infrared logarithms, it has only three; and so on. 



^We shall use henceforth (( O )) to indicate that the expectation value of an arbitrary 
operator O is taken in the leading logarithm approximation. 

^One can either view this kinetic contribution as part of the stress tensor, in which 
case it requires no composite operator renormalization but it also fails to contribute at 
leading logarithm order. Such divergent, subleading logarithms are automatically removed 
by renormalization as a consequence of the cancellation of overlapping divergences, but 
they can leave calculable, finite remainders. (For worked-out examples see [^.) Or one can 
view the coincident kinetic term as a composite operator in its own right, which of course 
requires composite operator renormalization to produce a finite result. For an explicit 
example of how even the leading logarithm contributions to such composite operators can 
be divergent see section 4 of [9]. 



2.3 Quantum Gravity 

Working in the leading logarithm approximation effects a dramatic simplifi- 
cation in quantum gravity. Full quantum gravity is haunted by the yet un- 
known - and experimentally untested - nature of its ultraviolet completion. 
In perturbation theory this shows up as an escalating series of divergences 
and counterterms whose arbitrary finite parts degrade - but don't entirely 
eliminate [26] - the theory's ability to make predictions. In contrast, the 
leading logarithm terms of quantum gravity are largely dominated by the in- 
frared sector, whose preservation is an essential correspondence limit of any 
ultraviolet completion. And the only counterterm that contributes at leading 
logarithm order is the 1-loop renormalization of the cosmological constant. 

This crucial insight concerning counterterms derives from the rule given 
above. The cosmological counterterm is: 

IottG 

We can express 6 A as dimensionless contributions from one loop (5Ai), two 
loops {SA2) and so forth: 

SA = H^{ Ml X K^H^ + 6A2 X {K^H^f + . . . } . (29) 

The measure factor has its standard expansion: 

v^ = a^{l + ^Kh-u + ^K\huf + ...] . (30) 

It follows that the cosmological counterterm has the general form: 

A£i ~ H^{ 6A1 + 6A2 X K^H'^ + . . . } X fi:"/i"a^ . (31) 

The 1-loop term has the same number of undifferentiated gravitons per cou- 
pling constant as the bare interaction (122|) . so it contributes at leading log- 
arithm order. However, one can see that the 2-loop term - the one propor- 
tional to 6A2 - is suppressed by a factor of the minuscule parameter k^H^. It 
contributes to the first subleading logarithm term, and is completely absent 
from the leading logarithm approximation. The higher loop contributions to 
6A are even more suppressed. 
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The higher counterterms of quantum gravity can be organized so that 
their graviton expansions always contain differentiated fields: 



A£o = 



6G 



AA 
AA 



a 



IQttG^ 
R-DA 



R-{D-2)A 

2 , 



-9 



= pCf"""'C, 



pafiu ' 



-9 



and so on. As with 6 A, we can express 6G in terms of dimensionless, 
contributions 6Gf: 



(32) 

(33) 
(34) 

-loop 



D-2 



(35) 



(36) 



6G = K^{ 6Gi X K^H^ + 6G2 X {k^H^ + . . . } . 

Hence, the Newtonian constant counterterm takes the form: 

AC2 ~ K,'^H^[6Gi+6G2Xi^'^H'^ + . . .}{dhdh+Khdhdh+. . .} , 

When compared with the bare interaction (1221) even the 1-loop contribution 
has an extra factor of n"^!!^ for each undifferentiated graviton. Therefore, the 
Newtonian constant counterterm makes no contribution at leading logarithm 
order, and the same holds for all higher counterterms. 

The reason for this amelioration of the ultraviolet problem is that un- 
differentiated fields effectively lose their ultraviolet modes at leading loga- 
rithm order. One can only reach leading logarithm order if every term in 
the free field expansion that can contribute a factor of ln(a) actually does 
so. These factors of ln(a) come only from the infrared sector, that is, from 
H < k < Ha{t) in the free field mode sum iQ p!^. For example, the expec- 
tation value of two coincident free fields is: 

2D{D-3) 



n 



Xijit,x)xkeit,x) 



2D{D-3) 
{D + l){D-2) 



n 



X 



{D + l){D-2) 
1 



^i{k^t)j 



2^- 



5i{k5t}j 



71 2 
1 



Ti 



D-1 



■ 5ii5. 



ijVkl 



H 



u{t,k) , (37) 



D-1 



■ 6ii6. 



(47r) 



ijOke 



(UV) 



r(/^-i) 
r(f) 



2 ln(a) 



(38) 



To get the ultraviolet divergent constant (UV) requires the full range of 
integration and the complete form of the integrand. However, this divergent 



11 



constant is guaranteed to be of subleading logarithm order, so it holds no 
interest for us. To extract the finite infrared logarithm it suffices to retain 
only the first term in the long wavelength expansion of the mode function: 

We can furthermore restrict the range of integration to just the infrared: 



2^-2 ^^r 



D-1 „,D-1-^ 



rdkk^-'\\u{t,k)f -^ ii^4^ / " ^ , (40) 

Jh 27t^^ Jh k 

21n(a) . (41) 



27r^ 

The simplifications noted above do not apply to the correlators of dif- 
ferentiated free fields. Because these correlators cannot contribute infrared 
logarithms, the constants they contribute (times one power of a for every 
spatial derivative) come from all parts of the mode sum and from all terms 
in the mode function. It follows that the correlators of differentiated free 
fields must he evaluated exactly, with the dimensional regularization in place. 

It is well to close this section with a summary of facts about infrared 
logarithms in quantum gravity: 

• There can be at most one factor of ln(a) for each extra factor of the 
loop counting parameter GH^\ 

• The operator under study controls how many infrared logarithms occur 
at lowest order, and it may also cause the number expected at higher 
orders to be fewer than the maximum just noted; 

• The only counterterm that affects the leading logarithm order is the 
1-loop renormalization of the cosmological constant; and 

• Even at leading logarithm order the correlators of differentiated free 
fields must be computed exactly, with the (dimensional) regularization 
in place. 
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3 Effective Scale Factor Approximation 

The purpose of this section is to explain the Effective Scale Factor Approx- 
imation. We begin by simply defining the approximation. We then give a 
motivation for it. Finally, we comment on the potential problem of using 
expectation values rather than stochastic samples, and the closely related 
issue of spatial inhomogeneities. 

3.1 The Approximation 

The dynamical variable of synchronous gauge quantum gravity is the spatial 
metric, gij{t,x). On de Sitter background it is natural to express this as 
follows in terms of the graviton field hij{t, x): 

gijit.x) = a'^{t)exp Kh{t,x) , a{t) = e^* . (42) 

L Uj 

We have seen that the various components of the graviton field can be ex- 
panded in powers of the transverse-traceless free field Xij{'t,x), which is de- 
fined by expressions fl9l fT4l) . We also saw from fl38l) that infrared logarithms 
derive from the steady growth in the amplitude of Xij- 

The complete expansion of hij in powers of Xij reflects the full complex- 
ity of perturbative quantum gravity, including the ultraviolet. It is neither 
possible to retain all this complexity, nor would it even be desirable in the ab- 
sence of fully understanding the physical ultraviolet completion of the theory. 
What we seek instead is a simplified expansion that reproduces the leading in- 
frared logarithms. We propose that the Effective Scale Factor Approximation 
may provide such an expansion. The basic idea is that the graviton remains 
transverse-traceless and free, but propagates in the background geometry of 
an effective scale factor A{t) which is determined from the expectation value 
of the QoQ gravitational constraint equation of motion. More precisely, this 
approximation is defined by three statements: 

1. The full metric is the square of a (D-number, effective scale factor A{t) 
times the exponential of a transverse-traceless mode sum Hij[A]{t,x): 

g,,{t,x) — > A\t) exp[t,H[A]{t,x)l. . (43) 

2. The transverse-traceless mode sum Hij[A][t,x) is the same as the free 
field Xijit,^) with the de Sitter mode functions fITll) replaced by the 
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functions M[y4](t, k) that pertain for scale factor A{t) 



Hij[A]{t,x) = 



3. The effective scale factor A{t) is determined from the leading logarithm 
expectation value of the goo equation of motion ([3]), including the 1-loop 
renormalization of the cosmological constant: 



n 



Roo + ^R-^{D-2)[a+6A,k^H' 



n 



. 



(45) 



Note that the expectation value of any operator in the presence of a 
homogeneous and isotropic state such as the vacuum must itself be 
homogeneous and isotropic, even if the operator is not, so ( l45l) can 
indeed be regarded as an equation for A{t). 



3.2 Motivation 

In terms of the graviton's general decomposition i^ it will be seen that the 
Effective Scale Factor Approximation amounts to the following assumptions 
about the behavior at leading log order: 



hj{t,x) 
h{t,x) 



Hij[A]{t,x 
, 

[a2(t)_ 
(D-l) In 



mt)_ 



(46) 
(47) 

(48) 
(49) 



Because we do not yet know how to derive these assumptions from the leading 
logarithm approximation it is best to regard them, at this stage, as being in 
addition to it. However, we susepct that a proper derivation exists and we 
advance the following reasons for this belief: 



A key distinction between gravity and the various scalar models 

that have been studied is that a constant value of the graviton field has 
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no gauge invariant significance. Unlike scalar theories, gravitation can- 
not develop an effective potential other than the cosniological constant 
itself. There are certainly quantum corrections, but they must always 
involve derivatives, which reduces the number of infrared logrithms. 
This suggests that the leading logarithm theory may be "free" in some 
variable. A plausible candidate for this variable is hjj^{t,x) because 
it is the only component of the graviton that involves Xiji^i^) ^^ first 
order. 

By rotational covariance the transverse field hj{t,x) must involve a 
spatial derivative, which should preclude this component of the graviton 
from ever contributing at leading logarithm order. 

Note that constraint equations determine h^(t,x) and h{t,x) as func- 
tionals of hj^ . The full solutions are certainly not spatially homo- 
geneous. However, every spatial inhomogeneity means that a pair of 
transverse-traceless free fields has failed to contribute an infrared log- 
arithm, which is spatially homogeneous. Hence the leading logarithm 
result should take the form (HHHHD we have assumed. 



• The Qqq constraint equation of motion would ordinarily determine a 
linear combination of h^{t,x) and h{t,x), so it is natural to use it to 
infer the effective scale factor A{t). Taking the expectation value at 
leading logarithm order ensures consistency of the method. 

It will be seen that the Effective Scale Factor Approximation provides the 
possibility for a self-consistent test of the "null hypothesis" that infrared log- 
arithms make no significant corrections to inflationary quantum gravity [T9] . 
For if the null hypothesis were correct then either the effective scale factor 
should remain at its de Sitter value — and the graviton should remain free — 
or else there are corrections but they drop out of a gauge invariant measure 
of back-reaction. In sections 4 and 5 we show that neither supposition is 
correct. 

Even though the long wavelength limit of a graviton is gauge equivalent to 
zero, it will be seen that each term in the graviton expansion of equation ( [73^ 
contains two differentiated gravitons. So any effect is sourced by the nonzero 
values of these derivatives and would indeed vanish for a purely constant 
graviton. However, recall that the expectation values of differentiated free 
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fields cannot be infrared truncated, so one has instead: 

n\xijit,x)xkiit,x)\n) ~ H"" , (50) 

{n\^mX^At,x)^nXkeit,x)\n) ~ a''{t)H^ . (51) 

These are still only small corrections, but they get enhanced by infrared 
logarithms at higher orders and the final effect may become large at late 
times. 



3.3 VEVs versus Stochastic Samples 

Inflationary cosmologists sometimes object to taking expectation values to 
follow evolution of a homogeneous, mean geometry, as we do in (H^ . They 
maintain that a long period of cosmological evolution presents a sort of 
Schrodinger Cat Paradox in which the wave function of the universe decoheres 
into many different portions, which are no longer in good quantum contact 
with one another, and in which the geometry is not even approximately spa- 
tially homogeneous on super- horizon scales [SB]- There is no question that 
this view is correct at some level because the vacuum expectation value of 
the stress-energy tensor is homogeneous and isotropic (because the state is) 
whereas we perceive (indeed we are) inhomogeneities and anisotropies. The 
key issues which have not been quantitatively addressed, and which bear on 
the validity of the present analysis, are: 

1. How unreliable are expectation values? and 

2. How much spatial variation should one expect? 

We suspect that the reliability of expectation values depends upon the 
operator under study. For operators which average to zero, such as the den- 
sity perturbation, the entire result arises from the decoherence effect, so one 
makes an enormous mistake in ignoring it. Other operators — for example, 
the square of a scalar fleld — acquire a signiflcant homogeneous expectation 
value upon which spatial variations are superimposed. Any quantum fluctu- 
ation drives this sort of operator positive, so one might happen to inhabit a 
special region of the universe in which there is little effect for a long time, 
but there will sooner or later be a large effect. The expectation value of such 
an operator can correctly reflect the long-term trend everywhere in space, 
even though it misses variations from one region to another. 
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The example of ^"^{t, x) is not specious because the free expectation value 
of this operator is one of the two ways infrared logarithms arise in perturba- 
tive computations. An even more likely candidate for the general reliabilty 
of expectation values is the gravitational interaction between the infrared 
gravitons which everyone agrees are produced during inflation. The self- 
gravitation of these particles must slow inflation, even if their distribution is 
not perfectly homogeneous. And this self-gravitation must grow with time 
as more and more infrared gravitons are ripped out of the vacuum by in- 
flation. As discussed in section 1, the fact that the total kinetic energy of 
infrared gravitons grows like the cube of the universe's radius means that 
the entire universe rather quickly falls inside its own Schwarzschild radius. 
The actual gravitational potential grows like In (a), rather than a?, because 
causality delays effects from a patch of spacetime until it passes within one's 
past light-cone p^ . Because gravity is a weak interaction, even at the scales 
of primordial inflation, the self-gravitation of inflationary gravitons cannot 
be significant until the past light-cone has become enormous. So a signifi- 
cant back-reaction from quantum gravity can only occur as the cumulative 
effect of very many small fluctuations within the past light-cone of the ob- 
servation point. Because the past light-cones of nearby points must largely 
overlap, there cannot be much spatial variation within a few Hubble radii. 
There might well be significant variation on very large scales, but even this 
should not change the long-term slowing trend everywhere. So we believe 
that quantifying this trend through the use of expectation values does not 
represent a serious error. 

We emphasize that the reliability of expectation values is a quantitative 
issue rather than a qualitative one. We do not maintain that the actual 
universe is exactly homogeneous and isotropic, as its quantum average is. 
Rather we argue that using expectation values is unlikely to produce signif- 
icant errors for the specific back-reaction effect we seek to study. Far from 
ignoring it, the inherently stochastic spatial variation of back- reaction plays 
an essential role in our model of inflation: it is how one can get scalar density 
perturbations of about the correct strength without the need for a fundamental 
scalar. 

To understand what we have in mind, suppose that the leading infrared 
logarithms grow to become nonperturbatively strong and that they sum up 
to produce a graceful exit to inflation. As noted above, there may be large 
spatial variations on vastly super-horizon scales, but the general trend should 
be correctly described by the average geometry. Because the slowing effect 
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derives from an enormous number of small fluctuations generated over the 
past light-cone of a very large (~ 10^^) number of e- foldings, there can be only 
small spatial variation, even over regions as large as the currently observable 
universe. (Because the ~ 60 e-foldings after the end of primordial inflation is 
very much smaller than the ~ 10^^ we envisage to end inflation.) But there 
is some spatial variation; what is its likely strength? Because the leading 
infrared logarithms are spatially homogeneous and isotropic, perturbations 
around them must be subdominant. In other words, inhomogeneities must be 
suppressed by at least one power of GH'^. Because the leading logarithms are 
assumed to give an order one effect (stopping inflation) the flrst subleading 
terms should induce a spectrum of density perturbations of order GH^. This 
is what one wants! The power spectra for scalar and tensor perturbations 
of wave number k can be approximately expressed in terms of the values of 
the Hubble parameter and deceleration parameter at the time tk when that 
mode experienced flrst horizon crossing 



^s(fc)~^^^ ' -PTik) r-^ GH\t,) . (52) 

There is no way to recognize a factor such as 1/[1 + q{tk)], and this level of 
analysis, but the factors of GH'^itk) come out exactly right. 



4 Computing the Effective Scale Factor 

The purpose of this section is to solve equation fj45|) for the effective scale 
factor A{t). Because differentiated flelds must be treated differently than 
undifferentiated ones, we begin by isolating the derivatives. We then explain 
how the various indices can be factored out to reduce any 2-point correlator to 
one involving two massless, minimally coupled scalars. Next we give a brief 
description of how the equation can be formulated for a nonperturbative, 
numerical evolution. We close by working out the 1- and 2-loop corrections. 

4.1 Isolating the Derivatives 

In synchronous gauge ((700 = — 1, 5'oj = 0) the purely temporal component of 
the Einstein tensor G^,^ takes the form: 

11 1 

Goo = Roo - ^9ooR = g g'^g''^ gtjgke - giugn ) + - ^^"^^i? , (53) 
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where ^^ ^^R is the Ricci scalar constructed from the (D — l)-diniensional 
spatial metric gij{t,x). We now extract the effective scale factor: 



so that: 



G, 



00 



g^j{t,x) = A^{t)g^j{t,x) 



l{D-l){D-2)f^ + l{D-2)^r^-g 



(54) 



V 



+ 8 g'^g'^i hjlki - giklji ) + ttttt ^""-'^R 



where the spatial Ricci scalar is: 



(D-Dj^ 






^iji,k ~ ^tki,j 
gikje — gijM 



+ g'^g'^^g""'" 



2^2 



mjk^ nli -■- mij^ nkl 



+ t^t^g'^'' 



. gik.m gjt-,n . gij,m gk£,n gik,£ gjm,n ' gik,£ gmn,i 



(55) 



(56) 



(57) 



At this point we are ready to substitute the Effective Scale Factor Ap- 
proximation form (1431) for the metric. The fact that (jij is the exponential of 
a traceless field means that: 

(i) we can drop all derivatives of the form g^^dcjij, 

(a) we can convert the second derivative terms to a total derivative - whose 
expectation value must vanish - plus products of first derivatives: 



g^'^g gikji - gij,k 

~ g ,ij + g''''T''g""' [ siik/ gjm,n + gik,m gnij — gik,m gji 

As a result, the equation for A{t) becomes: 

1 







n 



G, 



00 



n 



2 7174 



1 



P-2) 






--n 



+ 



2^2 



Q 



gijg g-mn^ _ g.^^^ g . ^ ^^ _^ g .^ ^^ g^^. 



g'^g^^gjkgu 



Q 



Q 



(58) 



(59) 



(60) 



It remains to act the derivatives on specific if jj [A] {t, x) fields and isolate 
the expectation values of the differentiated fields. For this purpose, it is 
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useful to recall a standard formula for the derivative of the exponential of a 



matrix Mj^: 



de^^ = e*^ i (9M + ^ [dM, M] + ^ [[dM, M] , M 

I ^ . o. 



+ ... 



(61) 



We can define the term on the right hand side as the differentiated matrix 
contracted into a 4-index object B[M]iji:f. 



9(e^^)^. = [e''')^^xB[M]kj,,xdM„ , 



(62) 



which implies: 



B[M] 



ijkl 



^ik^N + —, 



Hk^tj 



2! 



5ikMf^ - Mik6. 



Hk^viij 



ik^l'j 



+ 



1 r 
3! 



5ik{M^),, - 2MikMej + {M%kk, 



+ 



(63) 



The 4-index symbol B[nH]ijk(^ is particularly useful for expressing deriva- 
tives of cjij = exp[KH]ij because the differentiated metrics in fl60l) are con- 
tracted into inverse metrics: 



9'^g^^9jk9u = B[KH]ijkeB[KH]jipgX k^HmHj 



:^,iJT^k£-mn- - 

9 9 9 9ik,m9ji,n 

-p.ijy.k(.-mn- - 

9 9 9 9ik,m9ni,j 



(64) 

B[KH]ijMB[KH]jipqg'^'^ X K'^Hu,rnHpq^n , (65) 

B[K,H]nikeB[KH]mjpqg''^ X K Hke,mHpq^n ■ (66) 



It follows that the equation for the effective scale factor can be written as: 



;P-2) 



(D-1)^-A-5A,k'H'' 



Q 



B[KH]ijkeB[KH]jij 



jipq 



Q)) x(Q 



1^ Hkl Hpq 



Q 



+ {^\l B[^H],,,eB[^Hypqg^^ - ^B[KHUkiB[KH]m,pqg'' 



Q 



.^(a 



2tt tt 

"^ J^kt,m J^pq,n 



Q 



(67) 



4.2 Factoring Out the Indices 

Since the mode functions u[A]{t,k) depend upon k only through its mag- 
nitude, we can use rotational invariance to reduce expectation values of 
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Hij[A](t,x) to those of a fictitious scalar mode sum: 



^{t,x) = 



d^'-'k 



- e{k-H)[u[A]{t,k) e^^-^a(fc) + (cc) } 



(68) 



(277)^-1 

This is achieved by writing the expectation value as an integral of the form: 



d^~^k 
— — ^^ f{k) kiky 



(69) 



and then extracting the indices through the familiar replacements: 

1 



nj2 rvj 



ki kj kk ki I >■ 
fcj kj fcfc ki kjji kji I > 



D-1^'' ' 



iD + l)iD-l) 
1 



(D + 3)iD + 1){D - I 



^ij ^ki + ^ik ^jl + ^ie ^jk 
Sij ^ki Smn + " " ' 



(70) 
(71) 
(72) 



iFrora the mode sum (jUj) and definitions (II OH 111) we can evaluate the 
coincident 2-point functions of interest: 
(i) The undifferentiated 2-point coincident function, 



Q 



Hij[A]{t,x)HM[A]{t,x 
d^'-'k 



Q 



(27r)^-i 
2D{D-3) 



iD + l)iD-2) 

2D{D-3) 
{D + l){D-2) 



u[A]{t,k) 



SilkSfM — 



nj(fcn^),j 



Ji{kOe)j 



5i{k5i)j 



SijSke 
D-1 

^ij^ki 

D-l 



D-2 
d^'-'k 



(27r)^-i 



u[A]{t,k) 



(a) The "temporal derivatives" coincident 2-point function, 



Q 



Hij[A]{t,x)HM[A]{t,x) 



Q 



(73) 
(74) 
(75) 

(76) 



2D{D-3) 
{D + l){D-2) 



Si{kSe)j 



D-l 



X {f 
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(in) The "spatial derivatives" coincident 2-point function, 



Q 



dmHij[A\it,x)dnHki[A]it,x) 
2D 



n ) (77) 

= {D + 3){D + l){D-l)iD-2) [ -(^ + l)^^:'-^'=^^-" + ^^^MmSn)i 
+ 4:5ki5i(mSn)j + iD'^-5)Si(kSi)j5mn - 4(-D - l)(5i)(fc(5^)(„(5„)(j X (llVv^f 

4.3 Nonperturbative Formulation 

The fictitious scalar expectation values in expressions (17511771) can be writ- 
ten as one dimensional integrals involving the mode functions u[A]{t,k) for 
general A{t): 






2D-2^{^)Y(D^)Jh 



dkk 



D-2 



2D-2^i^)Y^D^)JH 



dkk 



D-2 



u[A]{t,k) 



u[A]{t,k) 



2D^2^{^)Y(D^)Jh 



dkk 



D 



u[A]{t,k) 



(78) 
(79) 
(80) 



Some work remains to be done converting the expression for u[A\{t,k) ob- 
tained in [271 ES] to a form for which these integrals are useful, but this seems 
straightforward. In the remainder of this section we explain how (178ll can be 
used to compute leading log expectation values of nonlinear functions of the 
undifferentiated free field Hij[A\{t^x) such as those in the nonperturbative 
evolution equation ^7\\ . 

Let us begin with the expectation value of some analytic function F{x) 
of the fictitious scalar ip{t,x). Because the function is analytic we have: 



F{^) = E 



F"(0)(^" 



n=0 



n\ 



N 



Because the fictitious scalar is free we have: 
Thus, the expectation value of -^(v^) is: 



oo Tp2n 



F{V)) = E 



n=0 



F2"(0) 
2"n! 



X {^r 



(81) 



(82) 



(83) 
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In view of the identity: 



^ / dy e-y'y^^ 



{2n 



we reach the final form: 



dye-y' F{y.j2^)) 



(84) 



(85) 



We seek now the expectation value of an analytic function of the matrix 
of free fields Hij [A] {t, x) . Since they are free, any such expectation value 
will break up into sums of products of expectation values of pairs of these 
fields, possibly with different indices. Expression (1751) allows us to reduce 
the expectation value of any Hij pair to a tensor constant times (v^^). There 
are many integral representations analogous to (185|) for correctly representing 
the combinatoric factors. Rather than integrating over a traceless dummy 
matrix, we choose the simpler option of integrating over a symmetric dummy 
variable Qij and then enforcing tracelessness directly: 



1 



(^^) 



TT 



[rfQi,]e-«'^-«>^FMQ 



D-l 



AD{D-3) 
\\ {D + l){D-2) 



r 



(86) 



There is no need to continue working in D dimensions since we are only 
interested in the ultraviolet finite, leading logarithm contributions from such 
expectation values. We can therefore take D = A: 



F{H, 



Vy 



^ l[dQ^j\ e-^-'^- f{\q,, - Ui, Qkk 



vr 



-,1^- 



(87) 



It is also very convenient to transform Q to a diagonal form V via a 3- 
dimensional rotation TZ: 

Tlij = (dij - diOj'j cos{6) + Eijk Ok sin(6') + 9i 6j , 



Q 



nvn' 



Q^ 



Ttih 'Due TZ 



M I'^ji 1 



(89) 



This is advantageous because the rotation resides only on free indices in com- 
plicated functions such as the inverse metric and the 4-index object B[M]ijk(, 
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of expression 



exp 
B 



-kQ + -Tt{kQ)I 

o 



T^i-rrJljn X exp —kV+-1i{kV)I 



nQ--Ti{KQ)I 
o 



ijk£ 



I'^iml'^jnl'^kpl'^f.q ^ -D 



kD--i:i{kV)i 



and drops entirely out of the trace of Qij in the exponentiah 

Q^J Qij = I?n + ^22 + ^Is = T^l + ^2 + T^l ■ 

The measure factor in fl87p takes the form: 



, (90) 
(91) 

mnpq 



(92) 



(93) 



where J is the Jacobian determinant of the 6x6 matrix of derivatives of the 
six Qi/s with respect to the six variables {9\Vj). Because J is independent 
of 6^ we can compute it using the small angle approximation: 



T^ij — Sij + eijkOk + 0{9 



Then the transformation (l89l) gives 



Qii = Vi + 0{9^) , 

-(V2-Vs)ei + 0(6^) 



Q23 — Q32 
Q31 = Q13 



The Jacobian determinant is therefore: 



J 



iVi-V2)iV2-Vs)iVs-Vi 



Integrating out the angular variables gives a constant A^: 



[dQi 



i ij W i_ 



(94) 



(95) 



(96) 



We can find A^ by transforming to average and relative variables: 

V={Vi+V2 + P3) ; Ai = {V^ - V2) , A2 = {V2 - Vs) , (98) 
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so that: 



Jd^V \{Vi-V2){V2-V^){V3-Vi)\ e-'-^^-^^-^' 



,-ivf+vl+vl) 



g-§(A2+AiA2+Ai) 



dVe-"" / rfAi / rfAa AiA2(Ai + A2 

3 J — OO J — OO 

- / rfAi / rfA2 AiA2(Ai + A2 

We then change to sum and difference variables: 

a = Ai + A2 , /? = Ai - A2 
and evaluate the trivial Gaussian integrals: 



g-|(A2+AiA2+A2) 



(99) 



(100) 



d^V 



-{vl+vl+vD 






_ia2_ig2 



3vr 

^ • (101) 



Equation (llOip implies that the constant A^ in ( 1971) is: 



N 



Svr 



(102) 



4.4 The First Loop 



It is a simple matter to evaluate the right hand side of (1671) at one loop - 
that is, at order k^. Only the zeroth order parts of the inverse metric and 
the 4-index symbol (|63|) can contribute: 



-mn 



jrmn 



Bijkel^H] — > SikSi 



(103) 



Taking account of the fact that Hij[A]{t,x) is transverse and using expres- 
sions (I76tl77p gives: 



1 



(D-2) 



(D-l)—-A-6\iK^H^ 
A^ 



-(n 



ri 11 j T 11 n 



I] J-^i] 



n) + 



8^2 



Q 



/l- 11 ii rnlJ-', 



ij.m-^-^ij.m 



n) + 0{K^ 



DiD-3)\/^_,.,\ , 3 / 2|,v7 112 



'^V) + ^('^^l|Vv^||- 



0(k' 



(104) 
(105) 
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Expectation values of the fictitious scalar depend upon the effective scale 
factor A{t), which is itself a series in k"^. However, at this order we require 
only the well known de Sitter limit [6]: 

where: 

gf^dx^dx" = -dt^ + a?{t)dx ■ dx . (107) 

Substituting in fllOSp gives: 
D-2 






k'^H^ T{D + l) (4-3D)(D-3) 



0{k^) . (108) 



(47r)f r(f + 1) 16 

Since all factors in the above equation are finite, we can also set D = A: 

^^^n,H^-5A^K^H^ = -^^^ + 0{K^) . (109) 

A^ Stt^ 

Furthermore, we express (11091) more simply as an equation for the square of 
A/ A, and we can more fully anticipate the form of the next correction: 

^ = H4l + ^5A,K'H^-^ + 0{K'HHn{a))\ . (110) 

At this point we should comment on various other 1-loop computations. 
Had we chosen 6Ai to null the 1-loop correction to {A/Ay the result would 
be: 

/synchronous StT^ 

A previous computation in a covariant gauge found [ID]: 

(6Ai] =-7^ . (112) 

V / covariant ZTT 

Finelli, Marozzi, Venturi and Vacca worked in synchronous gauge but they 
used an adiabatic regularization to get [H]: 

Ml) = ^ . (113) 

Vfmvv 19207r2 ^ ^ 
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These differences reffect the well-known fact that counterterms can depend 
upon the choice of gauge and upon the regularization technique. The physi- 
cally significant fact is the universal agreement that any 1-loop shift in [A/A)"^ 
can be absorbed into 6Ai. 

One can also estimate the result from just infrared gravitons. A very 
early computation by Ford got [42j: 



(M,) 



= +4 • (114) 

Ford 7r2 ^ ^ 



By instead modeling each graviton polarization as a massless, minimally 
coupled scalar and restricting to the infrared, we find [18j: 

V /infrared IGtT 

The reason (11 111) and (11151) differ in sign highlights an important aspect of 
the leading logarithm approximation. Expression (11051) reveals that the 1- 
loop effect consists of a sum of the expectation values of the squares of two 
operators. If one makes an infrared truncation of the resulting mode sums 
this must give a positive effect, which could only be nulled by a negative 
value of 5Ai. However, that is not the correct way to work at leading loga- 
rithm order. It is always the full quantum field theory which defines results; 
the restriction to infrared modes is only justified when extracting an infrared 
logarithm. The differentiated operators that appear in expression (I105P can- 
not produce infrared logarithms. What they give instead is constants, and 
the values of these constants derive from the ultraviolet as well as the in- 
frared. The correct procedure is to compute the expectation value from the 
full mode sum, with the regularization in place. When this is done one can 
see from expression (11061) that the expectation value of ||Vv?|P is actually 
negative, and large enough that it predominates over the positive expecta- 
tion value of 0^. Of course there is no mystery about the expectation value of 
a square giving a negative result; the automatic subtraction of dimensional 
regularization has subsumed the positive, power law divergence - which does 
not, in any case, contribute as vacuum energy - and left a finite, negative 
remainder. 

Whether or not we make the choice (II lip to null the 1-loop contribution 
to A/A, it is still constant. This allows us to evaluate the derivatives of 
expression (I105p to higher order using the de Sitter result (11061) with H 
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replaced by Aj A: 



2 -2 

K if 



3'^' 4i + OU') 



i^(-^iiv^r 



327r2 A4 



(116) 
(117) 



327r2 A4 

Since the 1-loop correction to A/ A is constant we can actually ignore it 
compared to the leading logarithm terms that derive from the undifferenti- 
ated factors of nHij[A\{t,x) in the inverse metric and in the 4-index symbol 

B[KH]ijki'- 



2 -2 



iil-'w^^r 



327r2 



0{k^ 



+ OU"" 



;il8) 
;il9) 



4.5 The Second Loop 

We are now ready to evaluate expression fl671) at 2-loop order. A useful 
identity is the expansion of two contracted 4-index symbols: 



Vt 



B[KH]ijMB[KH] 



jtpq 



Q 



Sqk Sip + 



K 
12 

Oak 0£v + TT 



n 



<^,fc(i^')£p - 2H,, H,p + iH'),k 6ep fi )) + 0(«:^) , (120) 
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4:Sak 5pr, — 5„P 6kn — San <^, 



Jqk Oep — Oqi Okp 



qpOke 



2 2 



OU' 



'1211 



Contracting (11211) into (1751) and then making use of (IllSp gives the first of 
the three expansions we require: 
1 
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(123) 



It is by now clear that we are doing an expansion in powers of a dimen- 
sionless time dependent parameter which may as well be named: 
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(124) 



It is also clear that every term in our expansion will involve one factor of xH. 
(Higher derivatives can and do appear in subdominant logarithm corrections, 
but they cannot occur at leading order because every time derivative elim- 
inates a factor of t = ln(a)/if . This is another way of seeing that only a 
renormalization of the cosmological constant is required at leading logarithm 
order.) With this notation the first of our expansions fll23p assumes the 
simple form: 



1 



Vt 



B[fiH]ijki B[kH 



jipq 



Vtl^iVL 



1^ Hki Hpg 



Q 



Hx[^ + ^x + Oix')] . (125) 



The inclusion of an inverse metric in relation (I12ip adds only a single 
extra term at the order we are working: 
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B[KH]ijke B[KH]jipq g"^ 
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Vt 
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~r ^ „ Omn 
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14:5nk. 5fr, — 5nP 6kn — S„„ 6, 



>qk 



q£ iJkp ~ iJqp ^kt 



Sqk ^Ip ^mn + — ^4:6 gk Sip — 6g£ 6kp — Sgp 6k£ Smn X + 0{x ) 



5™n((«:V')) + 0(«:4),(127) 

(128) 



The overall factor of 6mn in (I128P gives a simple result when contracted into 
the derivative term (!77|) : 



^(» 
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D-l 

We can take the D = A limit of the above equation because it is finite: 
1 



A^ 



n 



1^ dmHkedmHpg 
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—Hx 
10 



Sk{pSq)e — - ^ke Spq 



0{Hxx) . (130) 
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Multiplying the two factors (1128^ I130p gives the second of the three expan- 
sions we require for equation (1^ : 



- n 



B[KH]ijM B[KH]jipq g'' 



Hxl- 



27 117 



Q 



16 80 



x + Oix"^)] . (131) 



The two factors of the final term on the right hand side of (167|) cannot be 
so usefully expanded in isolation of one another. It is better to keep them 
together so that transversality can be exploited. The result is: 
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Si{kSi)j — 



X X 
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Si(-iSi>\k — 26ikS 



H(jOe)k 



ik^jt 



HxU} + ^x + 0{x^)\ . 



Hx + 0{Hxx^) , (133) 
(134) 



Substituting expansions (I125p . (I13ip and (11340 in our equation ( 1671) for 
A/ A gives: 



4- = H^ + l6AiK^H'-Hxl.l + lx + 0{x^)\ . (135) 

A^ 3 12 5 J 

Taking the square root and retaining only leading logarithm terms gives: 

H + -5Kik'^H^-x\- + -x + 0{x^)\ . 
6 14 5 J 



A 
1 



(136) 



Integrating reveals corrections to the effective scale factor as a series in powers 
of the parameter (I124p : 

ln[A(t)] = ln(a) + — (5AiX-ix-— x2 + 0(x3) , (137) 

where we remind ourselves that: 
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(138) 



One might think that relation fll37p reflects an unfortunate background 
dependence of our formalism. In the sense that the equation was derived us- 
ing perturbation theory it certainly does depend upon the zeroth order result. 
However, the larger issue is whether or not it represents the perturbative ex- 
pansion of a background independent result. For proper consideration of this 
matter it is crucial to distinguish dependence upon the initial state from de- 
pendence upon the background. The leading logarithm approximation, and 
our entire physical picture, only applies to an initial state which suffers a long 
period of inflation. This is a wide class but not a universal one. Many initial 
states are so heavily loaded with gravitons that they suffer gravitational col- 
lapse before beginning to inflate. We have nothing to say about such states; 
our focus is instead on states which are initial empty enough that inflation 
begins. The classical evolution of such a state is to locally approach the de 
Sitter geometry we used as the background in (jlj). In that case there is a per- 
fectly background independent way of interpreting the factors of In (a) = Ht 
in expression (11371) : H^ is the invariant acceleration measured by geodesic 
deviation and t is the invariant time from the point of observation to the 
initial value surface. 

5 Invariant Acceleration Observable 

Equation fll37p shows that the effective scale factor experiences secular slow- 
ing at 2-loop order. However, this is not enough to conclude that observers 
actually experience slowing, quite apart from the validity of the Effective 
Scale Factor Approximation. The problem concerns the noncommutativity 
of two operations: 

• Forming an invariant measure of expansion; and 

• Taking the expectation value. 

The variable A{t) can be regarded as the one third power of the local volume 
factor in synchronous gauge. Had the metric been classical, the logarithmic 
time derivative of A{t) would indeed give the expansion rate. But Unruh has 
shown that one cannot infer physics this way from the expectation value of 
the metric [13|. Instead of taking the expectation value and then forming 
an invariant, the correct procedure is form an invariant observable from the 
quantum metric and then take its expectation value. Models of scalar-driven 
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inflation which seemed to show secular slowing at 1-loop order from the 
expectation value of the metric [lU US] show no such 1-loop effect when 
examined with an invariant expansion operator [l6l W7\ . 

Finding an invariant measure of the expansion rate in pure quantum 
gravity is more difficult than for scalar-driven inflation because one lacks the 
preferred coordinate system in which the scalar is homogeneous. A reasonable 
proposal is based upon using the equation of geodesic deviation to measure 
the local acceleration between objects released from rest at some point x'^ 
^_lj. In the synchronous gauge we have been using this quantity takes the 
form: 

where A* is the spacelike separation of two initially parallel timelike geodesies. 
To make 7(2;) a full invariant one needs to average over the initial sepa- 
ration vector A* and also the position x^. A C- number such as A' cannot 
transform as a vector because only the fields of a quantum field theory trans- 
form. We therefore make use of an old trick [48j to convert it into a local 
Lorentz vector using the vierbein field e\(x): 

A' — > e\{x)A'' . (140) 

It will be seen that this takes 7(x) to the form: 

l{x) = y^9^j - -^9^^9ikgjije\e^^n^n' , (141) 

where n is the C-number unit vector in the initial separation direction: 

A^ 
n'' = , . (142) 

Va^a^ 

We can now take the C-number average over directions using 

d^~''nn^rf = j- , (143) 

to attain a scalar quantity: 
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To achieve a full invariant we must multiply by J —g{x) and integrate. 
Because the universe was released in a prepared state at t = 0, we can weight 
this integral by an arbitrary function of the invariant time from the initial 
value surface. That time is not an operator at all in synchronous gauge, it is 
just the coordinate time t. Moreover, we may as well choose the weight to be 
a delta function selecting a particular time. That still leaves the integration 
over space at this time. However, because the initial state is homogeneous, 
as is the gauge, we can dispense with even this step, although we do still need 



to multiply by J—g{t,x). The same procedure has already been exploited 
in computing the 1-loop expectation value of an invariant 2-point function 
in flat space background [I8] . 

We now make the Effective Scale Factor Approximation laid out in rela- 
tions ( I43ll45l) 11 We first evaluate 7 with the substitution gij = A^gtj, then 
take account of simplifications resulting from the tracelessness of Hif. 
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A^^ dtA 



A g'^gjj 

A D-1 ' 2{D-1) 

4(D-1) 
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(145) 
(146) 



Let us now take note of the highly significant fact that the measure factor is 
a (D-number in the Effective Scale Factor Approximation: 



^-9it,x) -^ [Ait)f . (147) 

We can therefore dispense with it altogether in the expectation value: 
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A^ d^A 

A^^diA^A{D-l) 



n 



B[KH]ijke B[K,H]jij 
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X IVL 



1^ HkeHpg 



(148) 



where we have also used relation 



^Note that this entails the potentially invalid step of treating a part of the average 
geometry — the effective scale factor A{t) — as though it is the quantum geometry. The 
fully correct procedure is to form an invariant operator first and then take its VEV. Our 
use of the Effective Scale Factor Approximation will only be justified if one can derive the 
approximation at leading logarithm order. We feel this can be done but we acknowledge 
that it has not been done yet. 
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At this stage we substitute equation (1^ for the effective scale factor in 
dHHD: 
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We have encountered each of the three terms on the right hand side (with 
shghtly different numerical coefficients) in the equation for the effective scale 
factor; taking D = i and substituting expansions (I125p . (I13ip and (11341) 
gives: 
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H' + ^^ + 1-6A,k'H^-^Hx{i+x + 0{x')} . (150) 



IFtoiu. expression (I136p we see that the time derivative of A/ A is always 
subdominant, and we should choose 6Ai so that initially the observable equals 
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Our final result is: 
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7 
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HU1--('^]' Ha) + [n^H^ In^ 



8 V 27r 
Hence, the observable 7 measures a slowdown of the expansion rate. 



(151) 



(152) 



6 Epilogue 

During inflation quantum gravitation loop corrections are enhanced by fac- 
tors of the number of e-foldings since inflation began. These enhancement 
factors are known as infrared logarithms and there can be at most one such 
factor for each extra power of GH^ in perturbation theory. The set of terms 
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for which this maximum number is attained is known as the leading logarithm 
approximation. 

Quantum gravity is vastly better behaved in the leading logarithm ap- 
proximation than in general. The reason for this is that, to reach leading 
logarithm order, every undifferentiated free field must contribute to an in- 
frared logarithm, which precludes these fields from producing ultraviolet di- 
vergences. From the plethora of BPHZ counterterms which permeate full 
quantum gravity, only the 1-loop renormalization of the cosmological con- 
stant contributes at leading logarithm order. 

We have explored a dynamical assumption called the Effective Scale Fac- 
tor Approximation. Under this assumption the graviton remains transverse- 
traceless and free, but propagates in the background geometry of an effective 
scale factor A{t) which is determined from the expectation value of the goo 
equation. By extending the formalism to the first subdominant logarithm or- 
der one can even describe cosmological perturbations of the correct strength. 

Under the Effective Scale Factor Approximation all perturbative compu- 
tations are straightforward. We obtain explicit 1- and 2-loop results for the 
effective scale factor and for an invariant observable based on using the equa- 
tion of geodesic deviation to probe the local acceleration. It is even possible 
to give a nonperturbative, numerical formulation to quantum gravity within 
the context of the Effective Scale Factor Approximation. 

The Effective Scale Factor Approximation may not be entirely correct 
but it does illustrate how simple inflationary quantum gravity can become 
at leading logarithm order. It also serves as a serious test of the "null hy- 
pothesis" that there are no significant corrections from infrared gravitons. 
For if the hypothesis was correct, then either the effective scale factor A{t) 
would receive no corrections or else they would drop out when evaluating the 
expectation value of the invariant acceleration observable. However, neither 
of these possibilities occurs. 
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